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> ! 1 INTRODUCTION 



ABSTRACT 

We investigate possible systematic errors in the recent measur ement of the matter 
power spectrum from the Lyman-alpha forest by Croft et al. ( 2001 ). We find that 
for a large set of prior cosmological models the Croft et al. result holds quite well, 
with systematic errors being comparable to random ones, when a dependence of the 
recovered matter power spectrum on the cosmological parameters at z ~ 3 is taken 
into account. We find that peculiar velocities cause the flux power spectrum to be 
smoothed over about 100 — 300km/s, depending on scale. Consequently, the recovered 
matter power spectrum is a smoothed version of the underlying true power spectrum. 
Uncertainties in the recovered power spectrum are thus correlated over about 100 — 
300km/s. 

As a side effect, we find that residual fluctuations in the ionizing background, 
while having almost no effect on the recovered matter power spectrum, significantly 
bias estimates of the baryon density from the Lyman-alpha forest data. 

We therefore conclude that the Croft et al. result provides a powerful new con- 
straint on cosmological parameters and models of structure formation. 

Key words: cosmology: theory - cosmology: large-scale structure of universe - galax- 
ies: formation - galaxies: intergalactic medium 



The Lyman-alpha forest is perceived as being simple, and in 
this simplicity promises to become a key to many unsolved 
problems in cosmology. In the most rec ent d emonstration 
of the power of simplicity, Croft et al. (2001) presented a 



measurement of the matter linear power spectrum on scales 
which are nonlinear today - a measurement not easily re- 
producible by other means, such as galaxy surveys at low 
redshift. 

Not only was the phenomenon that they modeled sim- 
ple, but also the method that they used to recover the linear 
power spectrum was simp le - p erhaps excessively simple. 



2 METHOD 

Since our goal is to evaluate the accuracy of the Croft et al. 
( [200 1| ) method, we apply their method to a range of cosmo- 
logical models. We adopt the same linear transfer function, 
but we vary the amplitude and the tilt, and we experiment 
with different geometries and different Hubble constants. In 



For example, Croft et al. (2001) considered only one prior 
cosmological model in deriving the correction function that 
translates the observed flux power spectrum into a linear 
power spectrum. 

In this paper our purpose is to check whether their 
method is reliable by sampling over a range of prior cos- 
mological models. We also attempt to estimate systematic 
errors and their covariance. We begin in § 2 by m aking sure 
that we are able to reproduce Croft et al. (2001) results. In 
§3 we give our results. We close in §4 with an optimistic 
conclusion. 



2 we also report the effect of varying the linear power 
spectrum over narrow intervals of wavenumber. We use a 
standard Particle-Mesh method to simulate th e distr ibution 
the matter at z = 2.72. Following Croft et al. (2001) we as- 
sume that the underlying baryon density follows that of the 
dark matter. At small scales the baryon density is smoothed 
out compared to the dark matter density, but the scale at 
which this occurs, the "filtering" sc ale, i s much smaller than 



the range of scales of interest here flSnedin fc Hui 1998[ ) (we 



discuss this in detail in § 3.3.2). All our simulations have a 
box size of 20/i -1 Mpc (in comoving coordinates) and 256 
particles on a mesh of the same size. The gravitational force 
is calculated with a Green function that includes the Opti- 
mal Antialising filter (Ferrell & Bertschinger 1991), by sum- 
ming over three Brilloun zones. 

In calculating the synthetic Lyma n-alp ha absorption 
spectra we follow the the Croft et al. (2001) procedure in 
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Fig ure 1 . The flux power spectrum predicted by the Croft et 
al. ( [200 1[ ) fiducial model, c ompar ed to the observed flux power 
spectrum from Croft et al. (2001). The dotted line is the model 



averaged over 3 best realizations (roughly equivalent to 100 ran- 
dom realizations) while the solid line is averaged over 12 best 
realizations (roughly equivalent to 400 random realizations). 



precise detail: we compute the synthetic spectra for 1000 
lines of sight which cross the computational box at random 
angles (to avoid bias introduced for lines of sight parallel to 
the box sides) for an assumed value for the spatial homoge- 
neous photoionization rate F. We then compute the mean 
opacity f(r) and find the value of F that gives the desired 
val ue for the mean opacity (r = 0.349 except as described 
in 



3.1.3). We then compute the ID flux power spectra of 



the normalized flux for each line of sight using a Lomb pe- 
riodogram and average the results over 1000 lines of sight. 
We have verified that this number - 1000 lines of sight - is 
fully sufficient to obtain numerical results with the precision 
well below the observational error. 

Befo re we can test for systematic errors in the Croft et 
al. ( 2001) method, we must first check that we can reproduce 
their results accurately. 

In order to reduce the uncertainty in the mean flux 
power spectrum, an average of several different random real- 
izations is usually used. Normally, this procedure requires a 
large number of realizations, because the variance decreases 
only as the square root of the number of realizations. In- 
stead, we adopt a slightly different approach. Out of all ran- 
dom realizations, we first choose "good" ones, i.e. those that 
give a measured linear power spectrum close to the input lin- 
ear power spectrum. We choose 3 "best" realizations out of 
100 random ones, and an average over those three recovers 
the input power spectrum at least as well as a plain average 
over 100 truly random realizations. 

Figure |l| shows the flux power spectrum Pf(&) (shown 
as A 2 (fe) = k 3 P(k)/2-K 2 ) averaged over 3 and over 12 best 
realizations (equivalent of 100 and 400 random realizations). 
The difference between the two curves provides an estimate 
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Figure 2. The flux power spectrum from the Croft et al. (2001) 
fiducial model (solid line) and our reproduction of their result 
(dashed line). 



of about 1/2 of the uncertainty in the mean flux power spec- 
trum (because averaging over 12 realizations gives half the 
variance of averaging over just 3 realizations). We notice 
that this difference is significantly smaller than the statisti- 
cal error-bars of the observational data, indicating that our 
flux power spectrum is calculated with sufficient precision. 

Hereafter we use three "good" realizations per normal- 
ization for each cosmological model. Thus, each model re- 
quires at least six simulations (two different normalizations 
above and below the best-fit value), and more if our ini- 
tial choice for the two normalizations does not bracket the 
best-fit value. 

As e vidence that our method reproduces the Croft et 
al. (2001) results, we show in Figure ^| their fiducial model 
(courtesy Rupert Croft) and the flux power spectrum from 
our simulation of exactly the same cosmological model. No- 
tice that there exist a difference on large scales, which is 
most likely due to the smaller size of our computational 
box, but that it is smaller than the observational errors and 
so is unimportant. 



3 RESULTS 

3.1 Systematic errors of the recovered linear 
power spectrum 

In this paper our main concern is the effect that an assump- 
tion of a specific cosmological model makes on recovering 
the linear power spectru m fro m the flux power spectrum. 
Specifically, Croft et al. (2001) assume that the flux power 



spectrum Pf(^) is proportional to the linear power spectrum 
Ph(k) at the same wavenumbcr 



P F (k)=b 2 (k)P L (k) 



(1) 
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the bias factor b(k) being independent of, or at least insen- 
sitive to, the power spectrum Pi(fc). 

If the "bias factor" b(k) were independent of the linear 
power spectrum, this would be a direct analogy to the biased 
galaxy linear power spectrum. However, this is not the case 
- the "bias factor" b(k) depends on the amplitude of the 
linear power spectrum Pl (k) , so that the relation ([!]) is not 
a linear relation. Correctly, equation (pf) should be written 



P F (k) =b 2 [k,P L ]P L (k), 



(2) 



where we use the square brackets to underline that b is actu- 
ally an operator acting on the linear power spectrum Pr,(fc). 
Thus, the recovered "observed" linear power spectrum, 



P°L hs (k) = 



P F hs (k) 
b 2 \k,P L ] 



(3) 



actually depends on the assumed linear power spectrum, and 
this dependence appears as a systematic error which is not 
i nclud ed in the estimate of the error bars by Croft et al. 
fl2O0 

For a particular choice of Pl that fits the observed flux 
power spectrum, relationship (JiJ) can be linearized. For small 
deviations APl about an assumed linear power spectrum Pl 
we obtain 



AP F (k) 



E 1 ' 2 



b 2 (k,k')AP L (k') 



(4) 



where the matrix b(k, k') can be obtained by differentiating 
equation (Q). Even in this linearized for m thi s equation is 
more complicated that the Croft et al. ( 2001 ) assumption 

(§)■ 

Throughout this paper, when we refer to the "recovered 
linear power spectrum", we always mean P£ bs (fc) obtained 
via equation (M). The sole purpose of this paper is to inves- 
tiga te ho w the dependence of b on Pl, ignored by Croft et 
al. (2001), affects the recovery procedure. 



3.1.1 Dependence on the Hubble constant or spectral 
curvature 



An important part of the Croft et al. (2001) procedure is 
the translation from the spatial units used in a simulation 
to velocity units of the observational data. This requires 
adopting a value for the Hubble constant at z = 2.72 which 
is used to scale spatia l scale s to velocity space. Specifically, 
following Croft et al. ( 2001 ), we adopt 

Wtot.i = Hz (n + v pcc j/H EdS>:j ) , (5) 

where n is the physical position of a pixel i in the synthetic 
spectrum, v to t,i is the velocity position of the same pixel, 
w P ec,i is the peculiar velocity of a given pixel in the simula- 
tion (in the EdS cosmology), Hz is the adopted value of the 
Hubble constant at z = 2.72, and -ffEdS.3 — 720ft km/s/Mpc 
is the Hubble constant in the EdS cosmology at z — 2.72. 
Thi s pro cedure is identical to the one adopted in Croft et 
al. (2001): we have applied our procedure to the simulation 



data kindly provided to us by Rupert Croft, and we obtained 
numerically indistinguishable results. 

It is important to emphasize here that this rescaling of 
the Hubble constant only affects the transformation from 
position to velocity space. The primordial power spectrum 
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Figure 3. (a) The flux power spectrum and (b) the r ecovered 
linear power spectrum p° bs in the Croft et al. (2001) fiducial 
model (EdS cosmology; solid line), the same cosmological model 
with the Hubble constant rescaled to a value appropriate to a flat 
cosmology with Q m .o = 0.4 (just like in Croft et al.; dotted line)), 
and for the flat cosmological model with Q m ,o = 0.4 (dashed line). 
Thin solid lines in this and all following figures show the assumed 
linear power spectrum Pl for the underlying cosmological model. 



of fluctuations as a function of position adopted in the sim- 
ulations is not changed, but rescaling of the velocity units 
does result in the linear power spectra being different in 
velocity space. 

It can be shown that for pure power-law power spec- 
tra all the dependence on the Hubble constant cancels out 
exactly. However, this is not true for a general power spec- 
trum with varying slope. In order to test this dependence 
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we have run a low density flat cosmological model in addi- 
tion to the fiducial model. Both models have identical initial 
conditions, and both run in the matter-dominated universe 
(f2 m> 3 = 1 for both models). They differ only by the val- 
ues of the Hubble constant H = -Ho-\/^m,o(l + z) 3 ' 2 , which 
are different by a factor y/OA = 0.6. Croft et al. (pOOll) did 



not use a Hubble constant appropriate to their cosmological 
model. Rather, they ran the simulation in the EdS cosmol- 
ogy, but adopted the Hubble constant from a flat cosmology 
with Q m ,o = 0.4. We also show the so rescaled model in Fig. 
| (which is the model also shown in Fig. ^). As can be seen 
from Figure ^[ the three recovered power spectra differ by 
less than the random error of the measurement, and we thus 
conclude that the dependence on the Hubble constant (or, 
equivalently, on the rate of change in the local slope - but 
see § 3.2) is insignificant. 



3.1.2 Assumed effective equation of state 



Croft et al. (2001) point out that uncertainty in the "effec- 
tive equation of state of the IGM", i.e. in the relationship 
between the temperature and the density of the photoion- 
ized gas, is one of the dominant systematic errors in recov- 
ering the linear po wer spectrum. This relation is commonly 
parameterized as (Hui fc Gnedin 199£) 



T^T (l + S) 



7-1 



(6) 



where To an d 7 a re parameters and S is the gas overdensity. 
Croft et al. ( fcOOlj ) adopt the values of To = 15,000 K and 7 = 
1.6. However, recen t measurements of the effective equation 
of state at z ~ 3 flRicotti et al. 200$ |5chaye et al. 2000 



McDonald et al. 2001) suggest somewhat different values for 
these parameters: T = 20,000 ± 3,000 K and 7 = 1.3 ± 0.2 
at z — 2.72 (all three measurements taken together and 
weighted by their uncertainties) . We have run the six models 
with values for To and 7 in the observed range, and the 
results are shown in Figure ^[ We find that the effect of the 
slope 7 on the recovered power spectrum is small, below 
the random uncertainties of the measurement, and so we 
do not attempt to model it, but regard it as part of the 
total systematic error quoted in Table [j]. The amplitude To 
however does affect the recovered linear power spectrum, 
which scales approximately inversely proportionally to To, 



(7) 



P L (k; To) = P L (k; 20,000 K) ( 20 '° 00K 

Rescaling the recovered power spectra this way reduces the 
rms difference between various models to about 1%. 

This dependence is easy to understand: for a given den- 
sity and velocity distribution and a given photoionization 
rate, higher temperature means broader absorption lines, 
which in turn implies a lower mean opacity. When the mean 
opacity is renormalized to the fiducial value 0.349, the pho- 
toionization rate is adjusted downward, and a given value 
of the flux now corresponds to a lower value of the density 
and thus the lower amplitude of the power spectrum (since 
readjustment of the photoionization rate scales all densities 
by the same factor) . 

Because the value To = 20,000 K is favored by observa- 
tions, we adopt this value as our fiducial value in the rest of 
this paper. 
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Figure 4. The flux power spectrum ( a) and the recovered linear 
power spectrum (b) in the Croft et al. ( [200 1| ) model with different 
assumed effective equations of state: To = 15,000 K, 7 = 1.6 (solid 
line, Croft et al. assumed values), To = 20,000 K, 7 = 1.2 (dotted 
line), T = 23,000 K, 7 = 1.2 (short-dashed line), T = 17,000 K, 
7 = 1.2 (long-dashed line), T = 20,000 K, 7 = 1.4 (dot - short- 
dashed line), and To = 20,000 K, 7 = 0.9 (dot - long-dashed line). 
Tq = 20,000 K, 7 = 1.6 (short-dashed - long-dashed line). 



3.1.3 Assumed mean optical depth 

Another possible systematic error discussed by Croft et al. 



( |200 ID is the assumed value for the mean optical depth 
- the value used to normalize synthetic absorption spec- 
tra. They assumed a v alue o f r = 0.349 based on Press, 
Rybicki, & Schneider jl993[ ) (hereafter PRS) value. The 
PRS measurement actually gives tprk = 0.349iKn!jJ . A 



more recent measurement by McDonald et al. (McDonald 
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Figure 5. The flux power spectrum (a) and the recovered linear 
power spectrum (6) in the Croft et al. model with different as- 
sumed values for the mean optical depth: r = 0.349 (solid line, 
Croft et al. assumed value), r = 0.260 (dotted line), t = 0.285 
(short- dashed line), r = 0.315 (long-dashed line), r = 0.400 (dot 
- short-dashed line). 



ct al. 2000) (hereafter MCD) gives a somewhat lower value, 
Tmcd = 0.285 ± 0.025. In order to test all the parameter 
range suggested by observations, we have run four mod- 
els with the mean optical depth different from the fiducial 
value: r = 0.400 (PRG+Ict), r = 0.315 (PRG-la and about 
MCD+lcr), r = 0.285 (MCD value), and r = 0.260 (MCD- 
ler). The flux power spectra and the recovered linear power 
spectra for these four cases and the fiducial case are given 
in Figure ^ As one can see, the recovered power spectra 
scale with r. We find a somewhat weaker dependence than 



that reported by Croft et al. ( [2001 ), mostly because we only 
consider a range of optical depth around the observed val- 
ues. The following simple power-law scaling takes out almost 
all dependence of the recovered linear power spectrum on r 
(the rms difference between the rescaled power spectra is 
less than 1%): 



P L {k;r) = P L (k;0 



349) ( 



0.349 



(8) 



The physical origin of this dependence is also clear: 
lower r means higher densities for the same value of the 
flux. 



3.1-4 Slope of the prior power spectrum 

Another possible source of systematic error is the slope of 
the prior linear power spectrum Pt,(k). We test four models, 
with n — 1.0, n = 0.95 (the Croft et al. model), n — 0.9, 
and n = 0.7. Figure ^| shows the flux and recovered linear 
variance for these models. The n — 0.7 model does not fit 
the flux power spectrum as well as other models, but the 
fit is still statistically acceptable (x 2 — 11 with 11 degrees 
of freedom). In all cases however, the recovered linear power 
spectrum agrees well with the fiducial model, indicating that 
the method is insensitive to quite significant variations in 
the slope of the prior power spectrum. In other words, the 
effective bias factor b(k) in equation (Q) appears similar for 
all the model spectra, an encouraging result. 



3.1.5 Dependence on the matter density 

Similarly, the dependence on the matter densit y para meter 
at z=2.72, O m , 3 , was not tested by Croft et al. ( |200 1| ). This 
may not be an issue after all, because for any reasonable 
flat model with a cosmological constant, £l m ,3 = 1 to better 
than 5%. Nevertheless, we have tested the matter density 
dependence. We have run three open models with n = 0.7 
and n m , = 0.45 (fi m , 3 = 0.75), n m ,o = 0.2 (fi m , 3 = 0.5), 
and fl m ,o =0.1 (f2 m ,3 = 0.3) respectively, and the results 



of these simulations are presented in Figure |^. One can see 
that the recovered linear power spectrum does depend on the 
assumed value of the density parameter at z — 2.72, even if 
all three models fit the observed flux power spectrum. 

The dependence on the density parameter cannot be 
rigorously derived analytically, but the form found empiri- 
cally can be understood approximately on the basis of argu- 
m ents from linear theory. The Gunn-Peterson optical depth 
t (Gunn & Peterson 1965) at every spatial point i s inversely 



prop ortional to the du/dx along the line of sight (Hui et al 



1997), where u is the total velocity (Hubble flow plus pecu- 
liar velocity) and x is the comoving distance. In the linear 
regime 



— ~ — = aH(l 
dx dx 



fS), 

where / = dln_D+/dlna » ft ' 6 



(9) 



. The first term in the paren- 
thesis is the Hubble flow, while the second one is due to 
peculiar velocities. In the linear regime the Hubble flow of 
course dominates, but at z = 3 most of the Lyman-alpha 
forest is mildly nonlinear, with = k 3 Ph{k) j '(27T 5 ! ) ~ 1, as 
can be seen from Fig. [^ 

If we assume that different k values are independent 
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Figure 6. The flux power spectra (a) and the recovered linear 
power spectra (b) in four cosmological models with the different 
slope of the primordial power spectrum, as shown by thick lines: 
n = 0.95 (solid line, Croft et al. fiducial model), the same model 
with n = 0.9 (long-dashed line), n = 1.0 (short-dashed line), and 
n = 0.7 (dotted line). Thin lines in panel (6) show the prior linear 
power spect ra for these models and symbols show the data from 
Croft et al. (Eoo3). 



(which is not quite true, but good enough for our hand- 
waving arguments), we would expect that the power spec- 
trum for the optical depth r is approximately inversely pro- 
portional to (1 + /Afe) 2 . In other words, we hypothesize that 
the recovered linear power spectrum depends on cosmologi- 
cal parameters in the following way: 



P° L hs (k) 



Pl ct (k) 



1 + A fc 
1 + fa A t 



(10) 



Figure 7. The flux power spectrum (a) and the recovered linear 
po wer sp ectrum (b) in four different cosmological models: Croft et 
al. (2001) fiducial model (solid line), an open model with f2 m ,o = 
0.45 and n = 0.7 (long-dashed line), an open model with f2 m ,o = 
0.2 and n = 0.7 (long-dashed line), and an open model with 
£lm,o = 0.1 and n = 0.7 (dotted line). 



where we will call P£ ct (fc) the "factorized" linear power spec- 
trum, fa is the values of / at z — 2.72, and 

1/2 

A fc = 1 " ) (11) 



pr(fc) should 



and thus depends on P£ bs (fc). In this form 
be nearly independent of the value of the matter density 
parameter at z — 2.72 fi m ,3. 

Figure ^ now gives P£ ct (k) for the four models shown in 
0. As one can see, the factorization ( |i"(i| ) does account for the 
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Figure 8. The factorized recovered linear power spectrum in four 
cosmological models from Fig. [?]. Line markings are the same. 
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Figure 9. The "best fit" recovered linear power spectrum (solid 
line) and its 1 sigma systematic error bars (dotted lines). 



main dependence of the recovered linear power spectrum on 
cosmological parameters, although some differences remain. 
It is however not convenient to use in practice because it is 
nonlinear and implicit for P£ bs (fc). But we can notice that 
for Q m ,3 substantially different from 1, is about unity 
for the range of scales where the observational error-bars 
are small. It turns out that a comparable factorization can 
be obtained if we adopt a fixed value = 1.4 in equation 



PZ hs (k) = Pl c \k) 



2.4 



1-4/3 



(12) 



This fit is linear, so it is easier to use in joint parameter 
estimation, while it still recovers the density dependence to 
about 2.2% precision. 



Table 1. The recovered factored linear power spectrum P£ ct 



(km/s)- 1 
1.99 x 10~ 3 
2.59 x 10~ 3 
3.37 x 1(T 3 

4.37 x 1(T 3 
5.68 x 1CT 3 

7.38 x 1CT 3 
9.58 x 1(T 3 
1.24 x 1(T 2 
1.62 x icr 2 
2.10 x 10~ 2 
2.72 x 10~ 2 
3.55 x 10~ 2 
4.61 x 10~ 2 



pfct 

(km/s) 3 
3.92 x 10 8 
3.06 x 10 8 
1.47 x 10 8 
3.85 x 10 7 
1.51 x 10 7 
2.08 x 10 7 
7.42 x 10 6 
4.23 x 10 6 
2.34 x 10 6 
1.25 x 10 6 
5.79 x 10 5 
3.15 x 10 5 
1.51 x 10 5 



Random 
err. (km/s) 3 
5.02 x 10 8 
1.80 x 10 8 
0.63 x 10 8 
2.59 x 10 7 
1.00 x 10 7 



0.41 x 10 7 
1.40 x 10 6 
0.67 x 10 6 
0.43 x 10 6 
0.18 x 10' 



,6 

0.54 x 10 5 
0.23 x 10 5 
0.09 x 10 5 



Systematic 
err. (km/s) 3 
0.90 X 10 s 
0.38 x 10 s 
0.53 x 10 s 
0.82 x 10 7 
0.44 x 10 7 
0.26 x 10 7 
1.46 x 10 6 
0.49 x 10 6 
0.18 x 10 6 
0.08 x 10 6 
0.33 x 10 5 
0.20 x 10 5 
0.10 x 10 5 



3.1.6 The systematic error of the recovered linear power 
spectrum 

We can now summarize our results. The recovered linear 
power spectrum can be factorized in the following way: 

let , 



Pl hs (k) 



Pr(k)QnQ T Qr, 



(13) 



where 



Qt 

Qt 



2.4 



1 + T4/3 

20,000 K/T , 



(0.349/t) 



(14) 



and P£ ct (fc) is independent of anything else (just numbers) 
and is shown in Figure |^ and in Table |l| together with its 
"systematic" error-bars. We put the word "systematic" in 
quotes because we considered only a subsample of all possi- 
ble cosmological models, and our sampling of the parameter 
space is by no means uniform. We thus suggest that errors 



quoted by us can be considered as an e stimate for the sys- 
tematic error in the Croft et al. (2001) measurement, but 



the precise value of such an error and its covariance matrix 

may depend very much on th e specifics of the pr i or. 

The observational data ( | Ricotti et al . 200C ; 3chaye et 



al. 2000; McDonald et al. 2001) sets the value of Qt at 



Qt = 1±0.15 



(15) 



and the error on Qt is uncorrelated with the systematic 
error on P£ ct (fc). 

With the mean optical depth the story is more compli- 
cated, because there is no sufficiently accurate measurement 
of it. The PRS data give Q t ,prs = l.OOt^os: while the MCD 
data give Q t ,mcd = 1.18 ± 0.07. Combining the two values 
in quadrature (because the data sets they used are indepen- 
dent), we find 



Q T = 1.11 ±0.05. 



(16) 
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Figure 10. (a) The assumed prior linear power spectra for the 
fiducial Croft et al. model with the power reduced at a particular 
value of the wavenumber. The power spectra are offset vertically 
by 0.3 dex for clarity, (b) Their respective best-fit flux power 
spectra. 



And for the most conservative estimate, we can include the 
whole range of quoted numbers: 



0.9 < Q T < 1.25. 



3.2 Band-Power Windows 



(17) 



Measurements of flux power sample power not at a single 
wavenumber, but rather over a finite band of wavenumbers. 
The band power windows b(k,k') 2 in equation (Q) can be 
extracted by differentiating the flux power spectrum Pp(fe) 
with respect to the linear power spectrum Pi(fe). The shape 



of the band-powers emerges most clearly if they are scaled 
with the fiducial powers, so we define scaled band-power 
windows by 



5 In P F (k) 
d\nP L (k') 



b 2 {k,k')P L (k)/P F (k') 



(18) 



In order to measure the effective band-power windows in 
the present case, we ran four addition al mo dels with the fidu- 
cial power spectrum of Croft et al. (2001) reduced by 20% 
in four different wave-bands, corresponding to four k values. 
Linear power spectra for these models are shown in Figure 
|ic|a, and their respective flux power spectra are shown in 
Fig. px|b. The difference between the different models is com- 
parable to the residual uncertainty in the mean flux power 
spectrum (Fig. |l|), and thus we are only able to measure 
the difference between the various models with the reduced 
power to an accuracy of about 30%. A more accurate de- 
termination of the band-power windows would require an 
implausibly large number of simulations. 

As can be seen, the abrupt reductions in the linear 
power spectrum lead not to abrupt reductions in the flux 
power spectrum, but rather to broad depressions. Figure 
pr] shows the inferred band-power windows themselves, the 
derivative of the flux power spectrum with respect to the 
input linear power spectrum at a given wavenumber, both 
for the case when the amplitudes of the input power spectra 
are kept fixed, and when the amplitudes are adjusted to fit 
the observed flux power spectrum. Evidently the effective 
band-power windows are quite broad. 

Fig. [ll] indicates that the band-power windows are 
broader at smaller scales. Figure ^attempts to quantify this 
feature by showing the same derivatives as Fig. |ll], except 
as a function of r — 2n/k rather than k. To bring out the 
similarity of shapes, the four band-power windows from Fig. 
[ll] have been shifted horizontally so that their center points 
ro = 27r/fco coincide (the centers of the band-power windows 
fco do not necessarily fall on k' , but are always within one 
bin in k space, from which we conclude that this discrepancy 
is most likely due to the finite fc-space sampling) , and verti- 
cally to the same maximum. The four band-power windows 
do not have exactly the same width in real space. Fitting to 
a Gaussian of width a, we find that a as a function of k' can 
be approximated as 



a(k') 



25 km/s 



■yjk' x lkm/s 



(19) 



In order to illustrate that it is indeed peculiar velocities 
that are responsible for the smoothing of the recovered linear 
power spectrum, we show in Fig. [l2] the derivative of the 
flux power spectrum with respect to the prior linear power 
spectrum at k' = 0.0096 (km/s) -1 (an equivalent of the long- 
dashed line) with all peculiar velocities set to zero. In this 
case the peak is much narrower, and its width is explained 
by finite sampling of the fc-space alone. 

It is important to underscore here that it is the power 
spe ctrum itself which is smoothed, and not the flux. Croft et 
al. (2001) did consider the effect of smoothing the Lyman- 



alpha absorption spectrum (Fig. 10 in their paper) , but the 
effect that we discuss here is different - it is a correlation 
of power between the neighboring k values rather than a 
reduction in power on small scales. 
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Figure 11. The derivative of the measured flux power spectrum 
with respect to the value of the input linear power spectrum at a 
given wavenumber for four wavenumbers, as labeled. The curves 
provide estimates of the effective band-power windows of the flux 
power spectrum. Filled circles show the values of k' . Two cases 
are shown: (a) the amplitudes of the input power spectra are fixed 
to ag(z = 2.72) = 0.23 and (b) the amplitudes of the input linear 
power spectra are adjusted to achieve the best fit to the observed 
flux power spectrum. 



This effect is easy to understand physically. Let us imag- 
ine that the matter power spectrum is a delta function, i.e. 
the power is only nonzero at a given value fco. In the ab- 
sence of peculiar velocities the three-dimensional flux power 
spectrum would be non-zero only in a narrow range of scales 
around fco (we cannot claim that this range is innnitesimally 
small because the relationship between the linear matter 
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Figure 12. The derivative of the measured flux power spectrum 
with respect to the value of the input linear power spectrum at 
a given wavenumber for four wavenumbers and for the fixed the 
amplitude of the prior power spectra shown as a function of r — ro, 
where r = 2n/k. This is essentially Fig. H] with a transformed x- 
axis. The bold solid line shows what the ong-dashed line becomes 
when peculiar velocities are set to zero. The curves are scaled 
vertically to have the same amplitude at maximum. The thin 
solid line shows a gaussian with a = 200 km/s, which is a fit to 
the short-dashed line. 



power spectrum and the Lyman-alpha flux power spectrum 
is nonlinear; however, the solid line in Fig. |l2| does demon- 
strate that this range is less than the fc-space sampling in 
Croft et al.). In the presence of peculiar velocities a given 
structure in physical space would appear shifted along the 
line of sight in redshift space (which is where the flux is 
measured) , thus smearing power over a range of scales com- 
parable to the ID velocity dispersion. 



3.3 Missing physics 

In the previous subsections we have considered only one pos- 
sible source of the systematic error, namely the choice of the 
prior linear power spectrum. Another possible uncertainty 
comes from the fact that no simulation includes all of the rel- 
evant physics. We can identify at least three major physical 
ingredients that are missed by pure PM simulations: 

(i) gas pressure; 

(ii) inhomogeneities in the ionizing background; 

(iii) shocks in the gas. 

We discuss these effects in that order. 



3.3.1 Gas pressure 

The role of gas pressure is well understood. In the linear 
regime it can simply be expressed by the filtering scale of 
baryonic perturbations, k,F'- 
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Pga s (fc) « PDM(k)e- 2k2/k i 



(20) 



where is a complica ted function of the entire thermal 
history of the universe ( Gnedin fc Hui 1998 ), and equation 
( pp| ) is a good approximation to the exact solution to the 
linear theory equations for k < 0.8kp. There are also strong 
reasons to believe that the filtering scale des cribes the effec t 
of gas pressure even in the nonlinear regime ( Gnedin 2000b ) . 

R ecent measurement s of the thermal history of the uni- 
verse fl Ricotti et al. 2000| ; |Schaye et al. 200C| ; |McDonaldei 
al. 200 j|) give kp — (35 it 5)h Mpc -1 , significantly smaller 
than the range of scales k £ 3/iMpc" 1 measured from the 
Lyman-alpha forest power spectrum (the reason being that 
pressure filtering, always lagging behind the growth of the 
Jeans mass, is subdominant to the thermal broadening of the 
absorption spectrum). Thus, the pressure filtering produces 
only about 2% correction to the recovered linear power spec- 
trum, and it is unmeasurable at the current time. This is, of 
cou rse, in a full agreement with the conclusions of Croft et 
al. (pOOl. 



This may explain larger numbers for the baryon densit y typ- 
ically found in such simulations (Haehnelt et al. 2001). 



One of t he poss ible concerns with this conclusion is that 
the Gnedin (2000a) simu lation was terminated at z = 4, 
whereas the Croft et al. (2001) measure the linear matter 
power spectrum at z — 3. 

However, observationally we know that the following 
three quantities do not change by more than a factor of two 
between z — 4 and z = 3: the mean photoionization rate 
( |Lu ct al. 1996| ; |Cooke et al. 1997|; |5cott et al. 2000| ), the 
star formation rate (Steidel et al. 199E), and the mean sep- 



aration between star forming galaxies (deduced from con- 
stancy of the mean UV luminosity density and the star for- 
mation rate). This immediately implies that the proximity 
regions of star forming galaxies occupy approximately the 
same fraction of the total volume at z — 4 and 2 = 3, and 
thus expected fluctuations in the ionizing background are 
comparable at z = 4 and z — 3. We thus can be confident 
that our conclusions about the role of fluctuations in the 
ionizing background at z — 4 also hold at z — 3. 



3.3.2 Inhomogeneities in the ionizing background 

Inhomogeneities in the ionizing background can come from 
two classes of sources: large scale inhomogeneities from 
quasars, and small scale inhomogenei ties from galaxies . The 
former is known to be unimportant (Croft et al. 1999), but 



the effect of the latter can only be estimated from advanced 
numerical simulations which include effects of t he rad iative 
transfer. Here we use a simulation from Gnedin ( p000a| ). We 
follow the Croft et al. (2001) procedure for the last output 
of this simulation at z — 4, but we analyze it two different 
ways. In the first case we use the neutral hydrogen fraction 
from the simulation data (which reflects the non-uniformity 
of the ionizing radiation background) to compute a synthetic 
Lyman-alpha forest at z = 4. In the second case we as- 
sume that the ionizing radiation is uniform and compute 
the neutral hydrogen fraction from ionization equilibrium. 
The value of the background is calculated by normalizing to 
the mean opacity in the simulation, so that in both cases 
the mean opacity is the same. In both cases we adopt the 
gas temperature as given in the simulation. After computing 
the flux power spectra in both cases, we find that they differ 
by less than 1%. 



The Gnedin (2000a) simulation suffers from its small 
box size (4/i -1 Mpc), and therefore a 1% difference is proba- 
bly an underestimate, but it appears unlikely that it un- 
derestimates the role of inhomogeneities in the ionizing 
background by, say, a factor of 10. We therefore conclude 
that inhomogeneities in the ionizing background do not bias 
the measurement of the matter power spectrum from the 
Lyman-alpha forest. 

Surprisingly, while the power spectra in the two cal- 
culations are very similar, the actual value of the volume 
averaged photoionization rate differs by 20%, being larger 
for the inhomogeneous ionizing background. Thus, simula- 
tions that assume homogeneous ionization and adjust the 
photoionization rate to fit the mean opacity of the forest 
actually underestimate the photoionization rate by at least 
20%, or, equivalently, overestimate the baryon density by 
at least 10% (and probably by a significantly larger factor) . 



3.3.3 Shocks and other mess 

Obviously, non-gravitational effects like shocks can signif- 
icantly affect the distribution of the L yman-alpha forest . 
Fortunately, we have observational data ( Ricotti et al. 200C; 
Schaye et al 2000t |McDonald et al. 2001 ) that suggest that 
non-gravitational effects are not significant in the forest. 
First, the IGM temperature at z ~ 3 is quite low, only 
20,000 K, and is entirely consistent with pure photoheat- 
ing of the gas. Second, the narrowness of the cut-off in 
the Doppler parameter vs column density distribution of 
the Lyman-alpha lines suggests that the photoionization- 
induced "effective equation of state" has a small scatter, 
again fully consistent with the pure photoheating of the gas. 
Any realistic shock heating would destroy this tight correla- 
tion. 

Thus, while we cannot rigorously prove that shocks and 
other non-gravitational effects can be neglected in modeling 
the Lyman-alpha forest, they would require a rather special 
fine-tuning which closely mimics the work of gravity and 
photoionization. 



4 CONCLUSIONS 



The two main drawbacks of the Croft et al. (2001) paper 



are (a) that it misses the dependence of the recovered lin- 
ear power spectrum on cosmological parameters (eq. Jl3]]) 
and (b) that it misses the fact that the effective band-power 
windows of the flux power spectrum are broad. The broad 
band-powers must induce broad correlations between esti- 
mates of the flux power, and consequently also of the re- 
covered matter power spectrum. But as the whole, we must 
acknowledge that the measurement that they present is re- 
markable - it currently offers the best measurement of the 
linear matter power spectrum over a range of scales that are 
nonlinear today. 

The effective band-power windows can be approxi- 
mated as Gaussians with a dispersion Ak — k 3 ^ 2 /(2ir) x 
25 (km/s) 1 ' 2 in Fourier space: 
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i r°° (fc j[ ) 

pREC (fe) = 1 / dfc 'pTRUE (fc , )e 2Afc2 ( (21) 

where Pf EC (fc) is the recovered and P™ UE (fc) is the true 
underlying linear power spectra respectively. The smoothing 
is caused by peculiar velocities. 

And while the resulting strong covariance between esti- 
mates of power at different wavenumbers will make careful 
comparison of the data and cosmological models more com- 
plicated than simple \ 2 fitting, this measurement will un- 
doubtedly play an important role in securing accurate values 
of the cosmological parameters. 

We also notice in passing that residual inhomogeneities 
in the ionizing background do not affect the measurement 
of the linear power spectrum, but may substantially bias 
estimates of the cosmic baryon density from the Lyman- 
alpha forest data. 
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